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ABSTRACT 


The microwave brightness toni|>orature measurements for Nimbus-5 Klec- 
t-ically Scanned Microwave Radiometer (KSMR), and Nimbus E Microwave 
Spectrometer (NEMS) are used to retrieve the atmospheric water va|x>r, liquid 
water and wind speed by a quasl-stnttstlcal retrieval technique. It is sho\ n 
that the ' rightness temperature can be utilized to yield these parameters under 
various weather conditions. Observations at 19.35 GHz, 22.235 GHz and 31.4 
GHz were input to the regression equations. The retrieved values of these 
parameters for portions of two Nimbus-5 orbits are presented. Then com- 
parison between the retrieved paramenia and the available observations on 

the total water vapor content and the surface wind speed are made. The esti- 

2 

mated errors for retrieval are approximately 0.15 g/cm for w ater vapor con- 

n 

tent, 0.5 mg/cnr for liquid water content and 0.0 m/sec for surface wind speed. 
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1. Introduction 

In this paper we will discuss techniques for measuring the water vapor and 
liquid water content of marine atmospheres and the wind speed at the oceans 
surface. Using microwave remote sensing the importance of water vapor and 
wind speed measurements arc well known and obvious. However, the liquid 
water content of the atmosphere has not lieen routinely measured and thus its 
relative importance is not firmly established. It is, nevertheless, available 
as a by-product of the remote sensing approach discussed in this paper and may 
well prove useful. 

Tills remote sensing approach involves the measurement of emitted micro- 
wave radiation at several wavelengths and sorting out the various contributions 
to the measured brightness temperature. In a recent paper by Wilheit and 
Fowler (1977). a linear approximation method is used to analyze the wind speed 
and cloud liquid water content from a set of microwave measurements collected 
by an aircraft platform. Their data set was taken under a rather narrow range 
of subarctic conditions which permitted use of a simpler interpretation technique 
than is used here. In order to describe the techniques involved, we will first 
review the relevant aspects of microwave radiative transfer. We will then de- 
velop a statistical technique for retrieval of geophysical parameters from the 
measurements. Using this statistical technique, a data set involving observa- 
tions over a typhoon it. the Pacific Ocean will then be analyzed. 



The specific sensors applied here are the Nimbus K Microwave Spectro- 
meter (NEMS) and the Electrically scanned microwave Radiometer (ESMR), 
both of which are in orbit aboard the Nimbue-5 satellite. P?talled descriptions 
of these instruments can Ik* found in Wllheit ( 1972* and Staelin et al (1972). For 
present purposes it is sufficient to note that the ESMR is a scanning instrument 
which measures the upwelling microwave radiation at 19. 35 GHz while scanning 
f>0 to either side of the spa ecraft and has a spatial resolution of 25 km at 
nadir. The NEMS instrumen* ’dews only at the nadir, and obtains observations 
at 22. 235. 31. 4 GHz and at three frequencies in the 50-70 GHz moleculav Oxygen 
band <n« of interest in this paper). It has a spatial resolution of 200 km. 


.1. Microwave Radiative Transfer in Atmosphere 

In the microwave portion of the electromagnetic spectrum and at tempera- 
tures typical of the Earth and its atmosphere, the Rayleigh- Jeans approximation 
for black body radiation works quite well. This enables radiative transfer for an 
isotropic, nonscattering, nonreflecting medium to be expressed as follows: 



- < I s - 


T o> 


U) 


where T |( represents the intensity of radiation at a given wavelength and in the x 
direction expressed as an equivalent black body temperature, called brightness 
temperature, a is the absorption coefficient and T„ is the thermodynamic tem- 
perature of the medium. If this is integrated for a uniform slab of thickness, 
a along the direction of propagation, x. the result becomes, 
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where T„ is the brightness temperature incident on the slab and T D 

i n ou ‘ 

represents the emerging brightness temperature. 

If the expression, i a, <s quite small or if I ,, - T B is small, the slab has a 
very small effect on the intensity of the radiation; whereas, if the slab is opaque 
t a - ) the radiation emerging out of the slab is characterized by the thermo- 

dynamic temperature of the slab independent of the intensity of the incident 
radiation. In this study, the atmosphere is divided into 100 layers each 0.2 km 
'hick. Each layer is characterized by the temperature, pressure and concentra- 
tion of const uents of the center of the slab. In a manner analogous to equation 
1 the brightness temperature. T B radiating from a surface is given by 

T b=<»-OT + cT o . (3) 

Inc 

T n la the brightness temperature incident on the surface, e is the emissivity 

Inc 

of the surface and T n is the thermodynamic temperature of the surface. The 
incident energy must also be averaged appropriately over all angles which 
can be reflected into the viewing direction. Since most natural surfaces are 
quite rough when considered on the scale of microwave wavelengths (0.1 - 30 cm), 
we will use the Lambertian approximation (Born and Wolf, 1975) for this angjlar 
distribution. 

J.w/2 

Tjj(t?) cos 6 sin 6 d 6 , (4) 

O 

where T H ( ) is the incident brightness as a function of the incidence angle f J. 
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For cnmiHitationul simplicity, since the weighting function. (Sin cos ) w ithin 
the integrand has a maximum at * 45 c , we will approximate T„ bv 

B lnc 

T„ T u (t»=>45°). (5) 

Inc 

Thus, to determine this brightness temperature which would be observed in a 
particular situation, we must lx? able to characterize the emuslvlty of the 
underlying surface and the absorption coefficient ot the atmosphere as a function 
of their physical and chemical parameters. 

The emissivity of terrestrial surfaces in the microwave region of the elec- 
tromagnetic spectrum varies from 0.1)5 for must solid surfaces to about 0.4 for 
the ocean surface. It is this low emissivity and relative homogeneity of the calm 
• wean surface which provides the possibility of sounding the water vapor and 
1 juid water content of the atmosphere. 

The emissivity of a smooth ocean surface can be calculated from the dielec- 
tric constant through the Fresnel relations (Jackson, 11)02). The dielectric con- 
stant data used here are from Lane and Saxton (1952). At a given frequency, v, 
the real part of the dielectric constant, < ' , may be expresse ! as. 

♦ (t , - „) ( 1 + 2 t 2 ), (G) 

and the imaginary part, «■ ", as, 

( - t>T (C g * 6 b ) (1+(, 2 T 2 )+o . , (7) 

where = 27rv. The - , , jt and r values at a number of temperatures and 
NaCl concentration are given by Saxton and Lane (1952). And conductivity value 
as a function of temperature from the International Critical Tables (1928). These 
can be summarized by the analytical expressions 

ta) = 4.9. (8) 
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r is the temperature in Kelvin and N the NaCl concentration in grant 
normal. The form used for the temperature dependence of r was suggested on 
theoretical grounds by S;L\ton and appears necessary in light of the large tem- 
perature dependence of this particular parameter. The simple polynomial fit 
appears adequate for the remainder of these expressions. We will approximate 
sea water as a O.tiN NaCl solution. 

When the sea surface is roughened by w ind the situation becomes more com- 
plicated. Calculations by Stogrvn <19<>7) indicate that there is little effect due to 
roughness for viewing at nadir. In horizontal polarization from nadir the rough- 
ness increases the emissivity ;ind the effect increases with viewing angle out to 
at least <>0°. In vertical polarization the roughness decreases the emissivity 
out to about 55° incidence angle, above which the emissivity increases w ith 
roughness and viewing angle. The observations of Hollingcr (1970) only partially 
confirm the theoretical expectations but both agree that there is little or no ef- 
fect due to roughness for viewing directly at nadir which is the case considered 
here. 

In addition to roughening the sea surface; the wind also creates foam which 
can also increase the emissivity of the surface. Nordberg et al (1971) found 
that for a frequency of 19.35 GHz the foam coverage increases the emissivity 
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ut nadir by 4.3 x 1 o ~ * foi each meter per second tin- wind speed exceeded 7 m/sec 
with no effect at w ind speeds below 7 in /sec. Measurements have been made by 
Van Melle et al. (11)7.')) at .?.G9 CHI/, on in artificially generated foam; found an 
effect comparable to that of Nordberg et al. < ID' ,. More recently Webster et al 
(1D7G) have measured the net effect of wind on the surface emisslvity at several 
frequencies from 1.4 to 37 GHz at view angles of 0° (nadir) and 38°. They found 
little frequency dependence above 10 GHz and a somew hat weaker effect than found 
by Nordbcrg et al. On this basis we w ill in this paper assume that for the fre- 
quencies of present interest (19.35, 22.235, and 31.4 GHz) that each meter per 
second that the w ind speed exceeds 7 in/s increases the surface emisslvity by 
3.2 x lo * 3 . The determination of the effect of wind speed on surface emisslvity 

as a function of frequencies, polarization and view angle remains an important 
research problem. 

In the Earth's atmosphere at frequencies between 1 and lot) GHz only three 
constituents need to be considered: water vapor, liquid water and molecular 
oxygen. All other constituents normally found in the atmosphere have either a 
sufficiently low abundance or weak microwave absorption that they are ignorable. 
We will therefore consider only these three constituents.. We will restrict this 
.reatment also to non-raining situations to simplify the treatment of liquid water 
n atmosphere. 

Water vapor has transitions at 22.235 GHz and 183.3 GHz in addition to many 
transitions at higher frequencies. Using the treatment given by Gaut (19<i8) the 
absorption coefficient can be expressed as: 


a = a 
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where kK refers lo the contribution due to 22.235 GHz ab*or|ition line and can 
he written as: 


KKS 

\\ here 
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is the line width of the 22.235 GHz absorption line. 

* water vapor density, in g/m* 

P » Pressure, in Torr 

T *» temperature, in Kelvins 

v = propagating frequenty, in GHz 

refers to the eontribution from all other lines and can be written as, 
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For liquid water, as already noted, only non-precipitating clouds are 
considered in this study. Since in such clouds the typical particle si/e is much 
smaller than the wavelength of interest, the Rayleigh approximation (Gunn and 
East. 1954) should be quite adequate to "epresent the scattering by cloud par- 
ticles less than 100 micron. The absorption coefficient can be written as: 


Cl oud 


(km -1 ) = 


0. 188 me 


(15) 


(«' ♦ 2 ) 2 + <e“)* 

where m = cloud liquid water content of atmosphere, in gm/m ’. 

The dielectric data r ' and t ) used are those of Lane and Saxton (1952) and 


discussed in the previous paragraph for the special case of pure water. The 
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combination of the factor of v in thi t* expression of «* anti < results in a fre- 
quencx dependence of approximately v* lor the absorption coefficient of clouds. 
Since clouds are often supercooled below the 205 K lower temperature limit of the 
l.ane and Saxton dielectric data on w hich an expression for t' and t are based, 
It is important that we have chosen an expression for the relaxation time, r w hich 
has some theoretical foundation. 

The absorption by the complex of molecular oxygen resonances near o.5cm 
can be approximated for frequencies below 45 Gllz b\ treating them as a line with 
a resonant term and a non- resonant term. The absorption coefficient can be 
written as: 


- 0 3 P Av Av 

T 2 (v - v 0 )* ♦ Av J v J + Av! 


(1G) 


where v, * GO GHz and v is the pressure broadening in GHz. 
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P > 250 tor r 


A v J 1 . 55 * 1 0- P (^)" “ ^31, 0.4^19- P> ) 250 > P 2 , 9 t or r (17, 


i / . _\ 0.85 

1 155 * 10-* 


P < 1 9 t o r r 


Calculations based on this approximation agree well with those based on full 
Meeks and Lilly (1963) expression as modified by Lenoir (1968) for frequencies 
less than 45 GHz. 


III. Quasi-Statistical Retrieval Technique 

As can be seen from the previous section, the brightness temperature which 
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would Ik 1 observed upwclling from the atmosphere depends on mam geophysical 

variables. However, several of the variables (temperature, humidity, and liquid 
water content) are, in themselves, eontinuous functions of height. Therefore, 
this problem has an infinite number of degrees of freedom. In the most general 
sense, a retrieval o! geophvsicul variables from a finite set ol observations 
could never Ik* performed. We must resort to approximations ami constrained 
sol lions if we are to derive useful information. 

One technique for retrieving geophysical parameters is the statistical tech- 
nique which was applied, quite successfully, by Waters et al (li*7.'>) to the problem 
of a|>proximating the atmospheric temperature profile from microwave brightness 
temperature measurements in the molecular Oxygen band ( ~ 5 mm). In this 

technique, they lx?gin with a statistical data base that consisted of a large set 

of actual tempe soundings. For each of these soundings, they calculate 

the expected upwolling brightness temperature at all frequencies at which meas- 
urements are to be made. The temperature at various heights are then regressed 
against the brightness temperature to determine the most probable (if the statis- 
tics are Gaussian) temperature from each level. 

If there is little information about a given level in the bi'ghtness temperature, 
then the retrieved temperature w ill tend toward the a priori average of the statis- 
tical data set and the residual of the regression w ill be little smaller than the 
variance of the a priori data set. It is to l>e noted that information about a given 
level can occur in a given measurement both directly through the equation of 
radiative transfer and also indirectly through correlations inherent to the 
atmosphere itself. 

This technique may be faulted in that it contains any bias of the a prior i data 
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set. When the brightness temperatures are strongly related in a direct sense to 
the geophysical parameters being i.Hrieved. this would be of little consequence, 
but where the coupling is through statistical correlations among the geophysical 
observables or when the retrieval offers only marginal improvement over the 
a priori statistics then the retrieved values will tend strongly toward the com- 
monplace situations. This may well give average statistics which suggest a very 
good retrieval technique but which gives poor results in the extraordinary situation 

which is often of the most interest. We, therefore, have modified this technique 
somewhat. We generate an artificial data set subtending the entire expected 
range in all the geophysical variables at issue and force any correlations among 
them to be zero or at least much smaller than expected in nature. As a by- 
product, this relieves us the necessity of providing a true statistical data base 
for parameters which have been measured sparsely or not at all. 

In the present case, we wish to interpret measurements made at frequencies 
of 19.35, 22.235, and 31.4 GHz. On the basis of the microwave properties of the 
atmosphere, we would expect the brightness temperature to be affected strongly 
by the wind speed at the sea surface and the net content of water vapor and liquid 
water in the atmosphere and rather weakly by sea rurface temperature, cloud 
temperature and details of the vertical distribution of water vapor. Thus vari- 
ables such as details of the temperature structure of the atmosphere would be 

expected to be quite significant. Thus for our statistical data base, we select 
' alues from handbook atmospheres (Air Force Cambridge Research Laboratory, 
196F) each with its associated water vapor content as specified in Table 1. In 
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Table 2, we specify the values of sea surface temperature, and surface wind 
speed and in Table :t we give the cloud models used. The statistical data base 
then consists of all possible combinations of these four parameters. This 
would provide zero correlation among the four parameters. However, we also 
specify that wherever the cloud liquid water content is not zero, the relative 
humiditv must In? 100 at that altitude. This gives a weak correlation between 
cloud liquid water content and net water vapor and provides some additional 
variability in the vertical distribution of water vapor. 'Thus, for the present 
case, at each of three frequencies we calculate brightness temperature for 
every possible combination, within our set of, sea surface wind, sea surface 
temperature, atmosphere model and cloud model. There are a total ol 1296 
data sets generated which contain the following information: sea surface wind; 
net water vapor content; net liquid water content; T (19. 35 GHz); T (22. 235 GHz); 
T b (31.4 GHz). Convention;:! linear predictions can then Ik? generated for each 
of the first three of these based on the last three while adding appropriate in- 
strument noise to the brightness temperature from a random numlier generator. 
While the linear approach works reasonably well, the dependence of the brightness 
temperature on the atmospheric water components is not linear. Rather the 
optical depth is more nearly linearly related to the water components. If. for 
the moment, we treat the atmosphere as a single layer with an optical depth 
T = T vapor + T ii l |.n<i :UK * 0 thermodynamic temperature t,, while the surface 
reflectance is It and has a temperature t^. We will also make the approximation 
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that t ( ~ t and we will ignore the cosmic background. Then the upwelling bright- 
ness T b is given by: 


T„ = ( 1 - K e _(,4 ‘ 2>r ) t , (18) 

The 1*2 occurs because of the approximation to a Lambertian surface discussed 
earlier. The logarithmic iorm of equation Is ean be written as 

Ln(t I -T,)«Ln(Rt,)-<l *>7>T v , por ♦ (1 ♦‘7)r Uqu|-# U») 

We then note that the function Ln(t, - T„) is linearlv related to each of the 
components of atmospheric water. Since the 22. 235 GH/. and 31. 4 GHz channels 
are primarily sensitive to water vapor and liquid water we will use these loga- 
rithmic functions rather than tha brightness temperature for these two frequencies 
and use the brightness temperature itself at 19.35 GH/.. litis discussion is not 
at all rigorous, but it serves to describe the gross behavior of the relationships 
and the functions so derived will mitigate the non-linear behavior of the problem. 
Since we are dealing with brightness temperature much colder than the thermo- 
dynamic temperature of the surface or atmosphere, we are not particularly 
sensitive to the value chosen for t,, thus for present purposes we will set 
t, ■ 280 K. We then solve for the three geophysical parameter wind speed (W), 
cloud liquid water content (I.) and water vapor (P) by conventional linear regres- 
sion so that the parameters are solved for in the form 
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where the values o r the a's and b's are given in table 4. In table 4 we also give 
the residual of each of the regressions which serves as an estimate the accuracy 
witl. which the retrieval can be made and the a priori standard deviation of each 
of the variables which specifies the variability in the data base. 

IV. Results 

To demonstrate the retrieval technique, two orbits of Nimbus 5 ESMR and 
NKMS data were chosen. In order to utilize boib »’SMR and NEMS data for the 
retrieval calculations, these data must be matched up geographically. Since the 
instantaneous field of view (1FOV) for the NKMS radiometers are quite different 
from that <>t KSMK, it is necessary to sum up all the 25 km x 25 km ESMR IKOV's 
which are contained within the NEWS' 1KOV that is approximately 2l)0 km x 200 km. 
The brightness temperatures measured at 0.90 cm, 1.35 cm and 1.55 cm radiom- 
eters for portions of orbits number 2s*.»l and 2t>is are shown in Figures 1 and 2 

plotted as a function of latitude. These matched data sets were then input into 
equation 20 to yield the estimates of the surface wind speed, cloud liquid water 
content and water vapor content in the atmosphere, i’he results ol these compu- 
tations are also shown in Figures 1 and 2 along with the corresponding ESMR 
images. By showing rain and land areas as dark shades, a spatial perspective, 
on the data is provided. The specific data sets used here are from Nimbus-5 
orbits 2 St) l and 29is which passed near the eye of 'typhoon Billie at about 1500 Z 
on July 14 and lti, 1073 respectively. Figures 3 and 4 give the corresponding 
surface analyses (each three hours earlier than the satellite pass) along with 
relevant observations of wind speed and atmospheric water vapor content. 
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Examination of the water vapor retrievals shows first that this varies smoothly 
with little structure. More specifically, on July 14, it can be seen that the re- 
trieved water vapor of about 4.5 cm at 30 'N agrees almost exactly w ith the ob- 
servations. Further, the slight decrease in water vapor along the track to about 
25°N is also suppoited by the observations. Beyond that, there are no supporting 
water vapor observations for this orbit. On July lo similar agreement iB ob- 
served between water vapor observations and retrievals. There are no direct 
observations of cloud liquid water content but the areas with liquid water re- 
trievals below 1 mg/cm correspond to clear areas as indicates by thermal 
infrared images from the same satellite. Also the fact that the bulk of the re- 
trieved data are very near to 0 mg/cm suggests that there is no consequential 
bias in this retrieval. Examination of the wind speed observations shows most 
observations in the 15-ilU kt (7.5-15 m/sec) range away from the storm center 
consistent with the retrievals and no observations in the high w ind portions of 
the typhoon. 

V. Conclusion 

\' e have developed a qunsi-statistical technique suitable for retrieval of 
geophysical parameters from microwave radiometric observations. We have 
derived specific equations for retrieving the wind speed at the oceans surface, 
and the atmospheric contents of both liquid and water vapor from simultaneous 
observations at 19. 35 GHz. 22. 235 GHz. and 31. 4 GHz and applied these equa- 
tions to data taken by tne Nimbus-5 spacecraft and by comparing the results 
with available observations shown that the results are promising. 
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Hu' estimate* of liquid water content in this study could l*e biased upward 
in the presenci of precipitation, because the present regression analysis is 
performed on the non-raining data ensemble. Further study in this area is 
required to resolve this problem. 
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Handbook Atmosphere 

1. 1'. S. Standard 

2. Tropical 

3. Subtropical summer 

4. Subtropical winter 
Midlatitude summer 

0. Midlatitude winter 

7. Subarctic summer 

8. Subarctic winter 

9. Arctic winter 


Table 1 

Water Vapor Content in cm J ) 
1.00 
3.88 
4.23 
2.12 
2.03 
0.87 
2.07 
0.42 
0.20 
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Table II 


Sea Surface Temperature Values 

<k) 

1. 273.0 

2. 283.0 

3. 293.0 

4. 303.0 


Surface Wind Speed Values 
(m/see) 

1 . 0 

2 . 10 

3. 20 

4. 30 
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Table 1U 


Cloud Top 
ikrn) 

Cloud Bottom 
<kni) 

Liquid water ck*nuit\ 
(K m J ) 

*> 

1 

0.01 

o 

mt 

1 

0.2 

8 

m 

7 

0.01 

8 

7 

0.2 

6 

1 

0.01 

6 

1 

0.2 

8 

6 

0.01 

8 

(i 

0.2 


No Cloud 
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Tab It* IV 


Regression 

Coefficients 

a t 

-100N. 

a 2 

1.831 

a 3 

37.92 

b n 

2.330 

b 12 

GO. 81 

b .3 

7G.68 

b 2. 

-0.0024 

b 22 

-0.014G 

b 23 

-0.2941 

b 31 

-0.0479 

b 32 

-8.699 

b 3J 

2.421 


A Priori 


Mean 

Standard Deviation 

Estimated 



Error 

2.35 

1.64 

0.15 

0.019 

0.030 

O.OOG5 


(U) Water vapor content (g/eni* 1 ) 
(b; Liquid water content (g/cm 2 ) 
(c) Wind 'peed (in /sec) 


15. 


11.2 


G.6 
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Figures 

Figure 1. The measured brightness tem|>erature» and the ESMK image of 

Nimbus 3 orbit no. jttOl, .lul> 14. 1973, The estlmalea t*f the water 
vapor content, the cloud liquid water content and the surface wind 
speed along the subsatellite track are prerented. 

Figure 2. The measured brightness temperatures and the KSMIt image of 

Nimbus 3 orbit no. 291». July lt» . 1973. The estimates of the water 
vapor content, the cloud liquid water content and the surface wind 
speed along the subsatellite track are presented. 

Figure 3. The surface weather chart corresponding to the studied area for 

1200/ July 14. 1973. The percentage of relative humidity and water 
vapor content (g cm-') for several weather stations are shown for 
references. 

Figure 4. The surface weather chart corre spending to the studied area for 

1200/, July 1(3, 1973. The percentage of relative humidity and water 
vapor content (g/cm J ) for several weather stations are shown for 


reference. 
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Table I. 

Table II. 

Table III. 
Table IV. 


Tablet* 

Handbook atmospheres used in this study and their corres|H»nding 
water va|Kir contents. 

The values of sea suriaee temperatures and surface wind speed used 
in this study. 

The cloud models and their liquid water density used in this study. 
The regression coefficients for equation (20) and the residual ui the 
regressions. 
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Figure 1. The measured brightness temperatures and the ESMR image of 
Nimbus 5 orbit no. 2891, July 14, 1973. The estimates of the water vapor 
content, the cloud liquid water content and the surface wind speed along the 
subsatellite track are presented. 
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Figure 2 . The measured brightness temperatures and the ESMH image of 
Nimbus orbit no. 2918, July 16, 1973. The estimates of the water vapor 
content, the cloud liquid water content and the surface wind speed along the 
subsatellite track are presented. 
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Figure 3. The surface weather chart corresponding to the studied area for 
1200Z July 14, 1S>73. The percentage of relative humidity and water vapor 
content (g/cnr ) for several weather stations are shown for references. 
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Figure -4. 'Phi* surface weather chart corresponding to the studied area for 
1200/ , July U>, 1973. The percent a no of relative humidity and water vapor 
content tg cm-) for several weather stations are shown for reference. 





